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We consider the interaction of a weak gravitational wave with electromagnetic fields 



o 
o . 

in a thin plasma on a Minkowski background spacetime using the 1+3 orthonormal 
frame formalism. Because gravitational and electromagnetic waves satisfy the same 
dispersion relation, electromagnetic waves can be effectively generated as a result of this 
[j^ ■ interaction. In the case of the interaction with a static magnetic field, the amplitude 

Qs, • of the electromagnetic waves depends on the size of the excitation region in which the 

magnetic field is contained. It is argued that due to the presence of a plasma this process 
, can also lead to the generation of higher harmonics of the original mode. Estimates are 

given for this effect in the case of a binary pulsar and a cold electron plasma. It is found 
that the emmited radiation will lie in the radio frequency band. We also speculate on 
the possible relevance of this process on situations in cosmology, in particular whether 
this could be used to constrain primordial magnetic fields. 
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1. Introduction 



It is well known that electromagnetic (EM) waves can scatter off gravitational fields. For 
instance, DeWitt & Brehme (1960) considered raditation damping in a gravitational field, and 
Stephani (1982) gives a general discussion on the subject. One can distinguish between two main 
areas: 

(a) The analysis of effects due to time- independent gravitational fields, such as Schwarzschild or 
Kerr spacetime, in an astrophysical context (Macdonald & Thorne 1982; Daniel & Tajima 1997), 
or 
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(b) scattering of EM fields off time- dependent gravitational fields (Cooperstock 1968; Zeldovich 
1973; Gerlach 1974; Denisov 1978; Grischuk & Polnarev 1980; Demianski 1985). 

The latter studies (b) have shown that there can be a coherent interaction between linear 
gravitational waves and EM fields, which would provide an effective means of transferring energy 
from gravitaional to EM degrees of freedom. It is for this reason that subsequent work focused on 
using this mechanism as a way of detecting gravitational waves (Lupanov 1967; Braginskii et al. 
1973; Grishchuk & Sazhin 1975). 

Motivated by the discussion given in the references above, we return to the time - dependent 
problem in order to further investigate a number of effects which arise when a gravitational wave 
interacts with an EM field, this time with a plasma present. In this letter we focus on a simple 
example, in order to demonstrate the basic effect, and in order to maximise clarity and make 
physical interpretations simple, we will use the 1+3 orthonormal frame (ONF) formalism (see, e.g., 
Ellis & van Elst (1998) and references therein). This has the advantage of picking out a physical 
observer and splitting spacetime relative to her. 

The paper is organized as follows: In section 2 we briefly review the ONF formalism and rewrite 
Maxwells equations using our ONF, introducing effective (gravity induced) charge and current den- 
sities. Next, in section 3, we specify the metric to that of a weak one-dimesional gravitational wave, 
and calculate the correpsonding (effective) current denstites. Section 4 then considers the excitation 
of an electromagnetic wave by a gravitational wave in a magnetized thin one - component plasma. 
We deduce that the resulting electromagnetic field amplitude is proportional to the gravitational 
wave amplitude, the magnitude of the static magnetic field and the size of the interaction region. In 
section 5 we consider the interaction region exterior to a binary pulsar. This consists of a vacuum 
regime (swept out by the Poynting flux) and a regime of interstellar matter. In order to simplify 
our analysis we consider only regions in which the magnetic field is on average aligned with the 
axis of rotation of the binary pulsar. Making an analog to laboratory experiments we argue that 
efficient harmonic generation of the EM waves may take place, and rough estimates suggest that 
these harmonics might - for favorable choices of parameters - be observable by the proposed As- 
tronomical Low Frequency Array (ALFA) (Jones et al. 1998). Finally, our results are summarized 
and discussed in section 6. 



2. Preliminaries 

Suppose an observer moves with 4-velocity u a (a = 0, ...,3), where u a u a = 1 (we put c = 1). 
This observer will measure the electric and magnetic fields (Ellis h van Elst 1998) 

E a = F ab u b , B a = \ e a b c F bc , (1) 

respectively, where F a b is the EM field tensor. Here e a bc = u d Vdabc, where r} a bcd is the totally skew 
4-dimensional volume element, and r/0123 = |det(g a b)| 1//2 (which will be just equal to 1 in our ONF 
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case). With these definitions, we can write the field tensor as 

F ab = u a E b _ u b E a + ^bc^ _ ^ 

We now introduce the ONF: e a , a = 0, ...,4, where eo = u, so that, with the above splitting 
(2) of the field tensor, Maxwell's equations V\ ) F ab = poj a , V[ a F bc ] = read 

V • E = p E + p m /e , (3a) 

V • B = p B , (3b) 
e (E)-V X B = -j E - n j m , (3c) 
e (B) + V X E = -j B , (3d) 

where the "effective" (gravity induced) charge densities and current densities are given by 

Pe = -rv^-e^rV^, (4) 

Pb = -r%5^ + ^r%£ 7 , (5) 

4 = [-(r%p-T a po )E^ + T%E a -e a ^(T /30 B 1 + r%B d )]e a , (6) 

j B = (-(rv-rv^ + rV^ + e^^o^ + r 5 ^)]^, (7) 



while p m and j m are the matter charge and current densities, respectively. Here Y% c are the Ricci 
rotation coefficients, and we have introduced the three- vector notation E = (E a ) = (E 1 ^ 2 ^ 3 ) 
etc., and V = (e\, e2, 63). The dot- and cross-products are defined in the usual Euclidian way. 

3. EM wave excitation by weak gravitational radiation 

To first order, the solution to Einstein's vacuum equations can be written as (using the trans- 
verse traceless gauge): 

ds 2 = -dt 2 + (1 + h)dx 2 + (1 - h)dy 2 + dz 2 , (8) 

where h is the deviation from flat space satisfying \h\ <C 1, and we consider only the +- polarisation 
of the gravitational wave. The deviation h is a solution to the wave equation Oh = 0, i.e. h = h(z—t) 
(Misner, Thorne, & Wheeler 1973). In order to clarify the physical intepretation and to simplify 
the equations, we introduce a contravariant ONF 

eo = d t , ei = (1 - \ h) d x , e 2 = (l + \ h) d y , e 3 = d z , (9) 

such that the metric components become 77^ = diag (—1, 1, 1, 1). Thus the local rest - space orthog- 
onal to a synchronus observer takes Euclidian form. 
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For the gravitational wave (8), the charge densities (4) and (5) are zero. The current densities 
(6) and (7) have components 

£ = -i(-E 1 + B 2 )h = g, (10a) 

-\{E 2 + B 1 )h = -j 1 B , (10b) 



where h = dh/dz = —dh/dt. 



4. 1 - D example 

Here we consider the effects of an incoming gravitational wave on a static magnetic field in 
an inhomogeneous plasma, restricting our attention to a 1-D case. Thus, we assume the presence 
of a external static magnetic field B = B(z)ei, and study gravitational wave excitation of small 
amplitude plasma waves propagating in the z-direction. We take the gravitational wave to be 
monochromatic, i.e. h = hexp[ik(z — t)], and we let the frequency of the excited wave coincide with 
the driver, that is, E a = E a (z) exp(— itot), where to = k, and similarly for all other quantities. From 
Maxwells equations (3), using |-B(z)|, \E(z)\ <C \B(z)\ and equations (10), we immediately obtain 7 

d 2 \ - - ( - 1 dB \ - 

— + k 2 \E 2 + \uj^f m = -lk 2 B--ik—\he W (ikz) = S(z), (11) 

u 2 E 3 + k^ Jm = 0. (12) 

In order to calculate the current, we assume the ions to be immobile - that is oj » uj p i, where u) p \ is 
the ion plasma frequency. Using our ONF the equation of motion for cold electrons can be written 

mu a v- n = q (v°E + v X B) , (13a) 

where v a = (v°,v) is the electron four - velocity, satisfying v a v a = —1, and u a is the observers four- 
velocity. Linearizing the equation of motion (13a) around the unperturbed plasma state (i.e. E = 0, 
v = 0) and dropping all cross terms containting the three- velocity and the rotation coefficients, 
we obtain 

-imujv = q^E + v X B^j . (13b) 
Solving equation (13b) for the electron velocity we obtain 
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UJ 2 



= l^-^E 2 , (14b) 



7 We linearize such that terms quadratic in barred quantities are dropped. 
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where u> c = qB/m is the electron cyclotron frequency. Using j m = qn^v where no is the unper- 
turbed electron number density at rest, and inserting equations (14) in (11) and (12), we obtain a 
driven wave equation for the extra - ordinary EM wave 

°" + k 2 - Ak 2 ) E 2 = S{z) , (15) 



9 

where the plasma induced wavenumber mismatch Ak is found to be 



.2^,2 ^ 



Ak 2 = -A, • (16) 

Here uj p = (nog 2 je^m) 1 ! 2 is the plasma frequency and uj^ = (to 2 + 10 If 12 is the upper hybrid 
frequency. Assuming that the static field is localized in the region < z < a, the exact solution to 
equation (15), when Ak = constant, 8 is 

E 2 (z) = exp [ifci oc 2:] / exp [— ifci oc z'] S(z')dz' 
Jo 

+ exp[-ik\ oc z] exp[iki oc z']S(z')dz' , (17) 

J a 

where the integration limits have been chosen so that there are no waves coming into the excitation 
region. Here k\ oc = (k 2 — Ak 2 ) 1 / 2 is the local wave number of the electromagnetic radiation. The 
solution (17) also holds approximately when both \Ak\ |fc| and \dAk/dz\ \kAk\. We note that 
for a large excitation region, ka 3> 1, EM waves propagating along the positive z-axis can obtain 
much higher amplitudes than those in the opposite direction, provided the mismatch is small, that 
is Ak k. In the next section we will consider the regime 

u 2 p /uj c «w<w p «w c , (18) 

in which case the condition AA; <§; A; is fulfilled. Focusing on this approximation, we may neglect 
the negative propagating wave all together, and write the solution in the outgoing region z > a as 

E 2 (z,t) = E out exp[i(/c out 2 - ut)} , (19) 

where 

ikh f a fi , , fiAk 2 

tJ ^ )exp b^ 



t], nl - -~ J B( z)vxp I ) d: (20a) 



and k out = h oc (z = a). 9 This shows that the effective size of the "gravitational EM wave transmit- 
ter" is determined by either the extension of the magnetic field, provided the direction of the static 



Specifically this result holds in a vacuum region, and in this case the effect has been used as an argument for 
employing EM fields as gravitational wave detectors (see e.g. Grischuk & Polnarev (1980) and references therein). 

9 If we insist on a sharply defined excitation region < z < a, k out typically becomes imaginary, reflecting the fact 
that waves with lu < uj p does not propagate in an unmagnetized one - component plasma. 



-6- 



magnetic field does not oscillate, or the mismatch distance L = 2k/7rAk 2 , whichever is smaller. This 
is a consequence of the fact that extra - ordinary mode and gravitational waves satisfy nearly the 
same dispersion relation in the regime of consideration, which tends to make linear wave interaction 
coherent over large distances. 

5. Applications to astrophysics 

It is an interesting question whether it is possible to detect the generated EM waves during 
realistic conditions. An immediate problem is that for most binary astrophysical sources, such 
as binary neutron stars (Rasio & Shapiro 1999), uj/2it < 10 3 Hz, whereas the plasma density for 
typical interstellar matter correspond to uj p /2it > 10 3 Hz. Although an external magnetic field may 
decrease the cut-off frequency significantly below lv p for EM waves propagating perpendicular to the 
magnetic field, we cannot rely on that effect for EM wave propagation over interstellar distances. 
Furthermore, observation from the ground are inaccessible for frequencises roughly below 10 MHz, 
due to the ionospheric cut - off. However, as we will demonstrate below, it is too early to conclude 
that gravitationally induced EM waves cannot be detected. 

Let us consider a binary system consisting of two identical pulsars, 10 with individual masses 
M Q (for a general review of the physical properties of binary pulsars, see Phinney Sz Kulkarni 
(1994)). This leads to a Schwarzschild radius Rs = 2GM & /c 2 « 3 km. Furthermore, we assume 
that the separation distance between the two neutron stars is d = 20i?s- Using a Newtonian 
approximation, we find that the gravitational wave frequency is uj ~ 2 x 10 3 rad/s. This corresponds 
to a situation fairly close to collapse. Using the quadrupole radiation formula (see, e.g., Stephani 
(1982)), assuming two Newtonian point masses, the expected time before merging is of the order 
of 1CT 2 s. 11 The unperturbed magnetic field is taken as B{r) = B SUT {r^ uh /r 3 , where the surface 
magnetic field -B S urf is of the order 10 6 — 10 10 T (where the upper limit refers to magnetars), and 
fpuis is the radius of the pulsar. Thus the fall -off of the magnetic field is that of a dipole. 

Next we want to estimate the maximum electric field from equation (20a). To get a magnetic 
field that is essentially static in the frame of reference where the binary system is rotating, we must 
exclude the region closest to the pulsar surface from the region of interaction. Assuming that the 
electromagnetic excitation effectively starts at a distance r = 60i?s from the center of mass of the 
system, using the estimates presented above, together with condition (18), we find from equation 



1 'One might be hesitant to apply our results derived for a one - component plasma, to the electron-positron plasma 
that exists in the vicinity of pulsars. For our case, however, the presence of the plasma is not very important in the 
EM excitation region, due to the strong magnetization which keeps down the mobility, but instead the plasma affects 
the EM wave propagation further away from the binary pulsar, as to be discussed below. For a discussion of EM 
excitations by gravitational waves in a multi - component plasma, see Ignat'ev (1997). 

11 During the last part of the inspiral the radiation is far from monochromatic. We do not believe that this will 
have any significant effect on the physics discussed. 
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(20a) that 

£max ~ 1.5 x 10 6 h iQt B SUT{ V/m , (20b) 

where we have denoted the gravitational wave amplitude at the beginning of the interaction region 
(r = 60i?s) by hint- Here we have used that the effective interaction distance is determined by 
a fall off in the interaction efficiency together with a spherical attenuation of the electromagnetic 
wave, rather than a decoherence of the modes as discussed in the analytical calculations. Using 12 
h int ~ 0.001 and B suri ~ 10 8 T, equation (20b) gives 

£ max - 50 MV/m , (20c) 

at a distance r max ~ 120i?s- I n this region, the induced electron velocity is kept non-relativistic - 
to a good approximation - by the external magnetic field, since the electrical and magnetic fields 
are perpendicular, and we have v ~ E max / B(r max ) »6x 10 5 m/s. Using an 1/r spherical decay of 
the EM field for r > 

r maxi the electron velocity can also be estimated at large distances. In a distant 
region where the magnetic field is assumed to have decreased such that uj c < u, the electron velocity 
is of the order v(r) ~ qE ma _ x r max / (mtur) , provided the ultrarelativistic regime is not reached. The 
radius r re i from the source where the electron quiver velocity is moderately relativistic (which we 
define to be v = 10 8 m/s) can then be solved for. We obtain 

r re i ~ 7.5 x 10 9 km ~ 500 au . (21) 

Since the magnetic field from the binary system is weak at this distance (i.e. the condition uj c < co is 
fulfilled), we deduce that the electron velocity must become relativistic at some distance r satisfying 
fmax < r < r re i - the precise value depend on extension of the magnetic field. 13 

When the electron velocities approaches the speed of light the excited EM wave becomes highly 
nonlinear, as the ponderomotive and relativistic nonlinearities in the momentum equation becomes 
comparable to the linear terms (Shukla, Yu & Tsintsadze 1986). Several nonlinear effects may 
take place, such as parametric excitation of other plasma waves (exciting lower frequencies) and 
harmonic generation of the original mode. A detailed study of the evolution of the EM waves is 
by far beyond the scope of our investigation. However, we note that the degree of inhomogeneity 
of the background medium may be crucial for the outcome. From now on we assume that the 
nonlinear regime is reached comparatively close to the binary source, where the plasma medium 
can be characterized as strongly imhomogeneous. This means that the background parameters uj c 
and Up are assumed to change significantly during a single wavelength of the EM wave. Strong 
inhomogeneities generally tend to supress parametric excitations (Rosenbluth, White & Liu 1973), 
and thus we expect harmonic generation to be the more important process in our case. To get a 



12 Using a Newtonian approximation with d = aRs, r = (3Rs, we find that h ~ (2q/3) 1 . 

13 It should be noted that the 1/r 3 behaviour of the magnetic field holds to the light cylinder, outside of which the 
magnetic field falls off as 1/r (Rees 1999). 
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feeling for the possible magnitude of this effect, we can compare with laser experiments (Carman, 
Forslund, & Kindel 1981; Carman, Rhodes, & Benjamin 1981), where the relativistic velocities 
induced in an inhomogeneous plasma leads to repeated harmonic generation by a factor of order 
50 of the original laser frequency, and these harmonics have a comparable amplitude to that of 
the original frequency (see Bezzerides, Jones, & Forslund (1982) and Grebogin, Triphati, & Chen 
(1983) for a theoretical explanation of this effect). For the parameters chosen in our example, the 
electron velocity can be as large as in the experiments presented by Carman, Forslund, & Kindel 
(1981) and Carman, Rhodes, &: Benjamin (1981). The inhomogeneity parameter 1/kL, where L is 
the inhomogeneity scale length, is also of comparable magnitude to the same experiments, and thus 
it is not unrealistic to assume that a significant part of the excited EM energy can be converted 
to higher harmonics with a harmonic number between 40 and 50 of the original frequency, which 
correspond to a EM wave frequency in the long wavelength radio regime, but well above the typical 
cut-off frequency of the interstellar plasma. Such EM waves can propagate over astrophysical 
distances and suffer no other significant damping than spherical attenuation. 

Detection of such low frequecy signals cannot be done on earth, due to absorption and reflection 
of the ionospheric plasma. Thus it is evident that no existing radio telescopes can measure the 
electromagnetic emissions discussed above, but well developed proposals like the Astronomical Low 
Frequency Array (ALFA) is expected to measure electromagnetic signals in the interval 30 kHz to 
30 Mhz (Jones et al. 1998). Assuming that 10% of the energy gets transferred to harmonic numbers 
with a high enough frequency for detection (roughly harmonic number n > 60), and estimating 
the minimum detection level (due to noise originating from the solar system) to 1000 Jy (Jones 
1999), the maximum distance of observation becomes ~ 3 x 10 5 ly if we assume that the radio 
signal is evenly distributed in the frequency range 30-40 kHz. Thus we deduce that the source 
considered above can be observed by a system such as ALFA provided it is located to our own 
galaxy. 14 The theoretical estimates of the event rates of collapsing binaries in our galaxy are of 
the order of 10~ 4 yr -1 (Tutukov & Yungelson 1993; Yamaoka, Shigeyama, & Nomoto 1993) (for a 
review, see Schutz (1996)). 



6. Summary and discussion 

In this letter we have investigated the interactions between weak gravitational waves and EM 
fields, using the 1+3 orthonormal frame formalism. In this way we were able to give a clear 
presentation of the result that an EM wave can be generated as a result of the interaction of a 
gravitational wave with a static isolated magnetic field. Furthermore, in the presence of a plasma 



14 Furthermore, we expect (at least for neutron stars with stiff equations of state (Rasio & Shapiro f999)) that at 
the final stage of the inspiral process, the frequency and amplitude of the gravitational radiation will increase, making 
the process more effective over a short time - span. Also a larger distance of observation can be achieved if we assume 
that the harmonic structure of the EM wave is not destroyed during its generation. 
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we argued that higher harmonics of the original EM wave may be generated, and an astrophysical 
example, consisting of a binary system, has been considered. 

Compact binaries, although not too common, are interesting sources of gravitational radiation, 
since the theoretical predictions of their spectral charachteristics are believed to be well understood. 
It is therefore an interesting question whether such charachteristics would be detectable through 
the radio wave emissions described above. In our case, the harmonic generation mechanism must 
be analyzed in more detail in order for useful spectral predictions to be made. If this was to be 
accomplished, it would shed light on the possible nature of the radio source, and thus would provide 
us with an indirect observational method for gravitational radiation. 

In addition to binary systems, there are of course other simultaneous sources of magnetic fields 
and gravitational radiation, e.g. quaking neutron stars. One of the more extreme events of this kind 
would be the quaking magnetar, with a magnetic field of the order of 10 10 T. The gravitational 
wave frequency of such systems may also be up to one order of magnitude higher than in the 
example considered above (at least for the gravitational wave tu-modes (Kokkotas, Apostolatos &; 
Andersson 1999)). However, the physics involved in generating the gravitational waves tends to 
be more complicated, and thus the corresponding estimates of the amplitudes and frequencies are 
more insecure (Andersson 1999; Andersson & Kokkotas 1998; Kokkotas, Apostolatos & Andersson 
1999). It may still be of interest for further studies, since the large magnetic fields of magnetars 
could compensate for the expected lower amplitude of the gravitational waves. 

It is of interest to consider possible effects due to strong gravitational field deviations from 
Minkowski spacetime. According to Daniel & Tajima (1997) all frequencies, e.g. the plasma- and 
cyclotron frequencies, are redshifted by the lapse function in a spherically symmetric spacetime. 
Since the important quantity in the above calculations are the quotient of certain frequencies, the 
effect will persist, if we assume that the present situation can be approximated as spherical symmet- 
ric. Of course, the frequency of the emitted radiation will be redshifted (w C mmision/w re ception = 
although z will be less than 2% for the interaction distance 60-Rs used in the above example. 

In light of these results we will discuss the possible cosmological implications in a forthcoming 
paper (Marklund, Dunsby, & Brodin 1999), using the 1+3 covariant approach and tetrad methods 
which in a clear unified way deals with such problems. For example, we speculate that this may 
lead to an imprint on the Cosmic Backround Radiation (CBR) spectrum through the interaction 
between primordial gravitational waves and magnetic fields of cosmological origin. Therefore CBR 
measurements could in principle be used to constrain cosmological magnetic fields. 

M. M. was supported by the Royal Swedish Academy of Sciences. P. K. S. D. was supported 
by the NRC (South Africa) and the URC (University of Cape Town). The authors would like to 
thank the referee for raising the question of core/crust quakes in magnetars. 
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